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Abstract 



We study QCD with two colors and quarks in the fundamental representation at finite 
baryon density in the limit of light quark masses. In this limit the free energy of this 
theory reduces to the free energy of a chiral Lagrangian which is based on the symmetries 
of the microscopic theory. In earlier work this Lagrangian was analyzed at the mean field 
level and a phase transition to a phase of condensed diquarks was found at a chemical 
potential of half the diquark mass (which is equal to the pion mass). In this article we 
analyze this theory at next-to-leading order in chiral perturbation theory. We show that 
the theory is renormalizable and calculate the next-to-leading order free energy in both 
phases of the theory. By deriving a Landau- Ginzburg theory for the order parameter 
we show that the finite one-loop contribution and the next-to-leading order terms in the 
chiral Lagrangian do not qualitatively change the phase transition. In particular, the 
critical chemical potential is equal to half the next-to-leading order pion mass, and the 
phase transition is second order. 
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1 Introduction 



QCD with two colors and quarks in the fundamental representation of the gauge group 
is in many respects similar to QCD with three colors and quarks in the fundamental 
representation. Both theories are confining and asymptotically free, both exhibit a spon- 
taneous breaking of chiral symmetry, and both are believed to be in a different phase 
at high temperature and density. However, the mechanism of the phase transition at 
nonzero baryon chemical potential has to be different. The nucleon in QCD with three 
colors is a fermion with mass much larger than the pion mass. At finite baryon density 
a phase transition may occur if the Fermi-surface becomes unstable. Instead, for two 
colors, the lightest baryon is a boson with the same mass as the usual pionic Goldstone 
bosons. It is a diquark state. For a quark-chemical potential equal to half the mass of 
the lightest baryon a phase transition to a Bose-Einstein condensate of diquarks takes 
place very much like the transition to a pion condensate for increasing pion chemical po- 
tential in ordinary QCD. This phase transition can be described completely in terms of 
a low-energy effective theory for the Goldstone modes associated with the spontaneous 
breaking of chiral symmetry [p], §• This low-energy effective theory essentially relies on 
the symmetries of the underlying microscopic theory. In earlier work ^ |2|] this theory was 
studied at the mean field level. A phase transition from the normal phase to a phase of 
condensed diquarks was found at a chemical potential equal to half the pion mass. This 
phase is a superfluid phase with massless Goldstone bosons. One of them results from 
the spontaneous breaking of the U(l) group associated with baryon number. 

For QCD with three colors and fundamental quarks, the fermion determinant at 
nonzero chemical potential has a complex phase making first principle Monte Carlo simu- 
lations impossible. On the other hand, for QCD with two colors the fermion determinant 
remains real at nonzero chemical potential and Monte Carlo simulations can be performed 
for an even number of flavors || f|. We are thus in a unique situation where both ana- 
lytical results and Monte-Carlo simulations are available in the nonperturbative domain 
of QCD at nonzero baryon density. The hope is that this ultimately will lead to a better 
understanding of dense quark matter in QCD with three colors and fundamental quarks. 
What already has been achieved is that the mean field results have been confirmed by 
lattice QCD simulations of several independent groups ||. 

The same analysis can be made in several other cases. They can be classified according 
to the Dyson index of the Dirac operator of the underlying microscopic theory || Q|. For 
QCD with two colors and Nf quarks in the fundamental representation the Dyson index 
is (3 — 1 with Goldstone manifold given by SU(2Nf)/ Sp(2Nf). For QCD with two or 
more colors and quarks in the adjoint representation the Dyson index is (3 = 4 and the 
Goldstone manifold is given by SU(2Nf)/0(2Nf). The third case with Dyson index 
(3 = 2 is QCD with three or more colors and quarks in the fundamental representation. 
Although chiral perturbation theory is irrelevant for full QCD at nonzero baryon chemical 
potential, the low-energy effective theory for (3 = 2 can be applied to phase-quenched QCD 
J7j, i.e. QCD at nonzero chemical potential with the fermion determinant replaced by its 
absolute value. Phase quenched QCD has been reinterpreted as QCD at finite isospin 
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density || and both theories are described by the same low-energy effective theory. This 
theory has been generalized to include chemical potentials for the different quark species 
II |10). The theories for /3 = 1, (3 — 2 and (3 — 4 display a similar phase diagram at 
finite baryon chemical potential. They can be studied within very similar effective field 
theories @. It is also possible to write down effective theories for the Goldstone sector 
of microscopic theories with an imaginary vector potential but without an involutive 
automorphism such as the axial symmetry of the Dirac operator. The main difference is 
the structure of the Goldstone manifold which is determined by the pattern of spontaneous 
symmetry breaking. Such theories enter in the context of disordered condensed matter 
systems [Jll, [l^] and the effective theory and its phase diagram are similar to that of 



QCD-like theories. They are also relevant to QCD in 3 dimensions at nonzero chemical 
potential [|13]. As is the case in chiral perturbation theory at zero chemical potential 
flU] , p~5| , |7|, [T7|, [18], [19|, ^(J these mean field results can also be obtained from a chiral 



Random Matrix Model |21] at nonzero chemical potential. 



Up to now the above effective theories have only been analyzed at the mean field 
level H], H, [F[ H H, pL0[| . In this work we perform a one- loop calculation for the case 
(3=1. Our first goal is to show that the theory is renormalizable in both phases, i.e. 
that for an arbitrary background field the one-loop divergences can be absorbed into the 
coupling constants of the next order Lagrangian. Our second goal is to investigate the 
one-loop corrections to the mean field results for the phase transition. Our expectation 
is that the critical chemical potential of the next-to-leading order theory is equal to half 
the renormalized pion mass at this order. Related questions such as the effect of one- 
loop corrections on the order of the phase transition and the critical exponents will be 
addressed as well. 



The paper is organized as follows. In section 2 we give a brief review of the low-energy 
effective theory for QCD with two colors and nonzero chemical potential. The next-to- 
leading order effective Lagrangian is introduced in section 3. The one-loop corrections 
of the leading order effective Lagrangian are calculated for an arbitrary background field 
and it is shown that the divergences can be absorbed into the coupling constants of 
the next-to-leading order effective Lagrangian. In section 4 we calculate the one-loop 
corrections to the free energy in the normal phase and the phase of condensed diquarks. 
A Landau- Ginzburg model for the phase transition is derived from the effective theory in 
section 5, and its parameters are calculated from chiral perturbation theory. In section 
6, we discuss the number density, the chiral condensate and the diquark condensate. 
Concluding remarks are made in section 7. Some of the technical details are worked out 
in two appendices. In Appendix A we calculate the next-to-leading order corrections to 
the pion mass and the pion decay constant. In Appendix B the one-loop integrals for the 
baryonic Goldstone modes in the diquark phase are evaluated. 
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2 Low-Energy Limit of QCD with Two Colors 



In this section we review the low-energy effective theory of QCD with two colors and 
fundamental quarks. More details can be found in the original literature [ 22| [l], [J. 



2.1 The Lagrangian of QCD with Two Colors 

The QCD partition function is given by the ensemble average determinant of the Dirac 
operator. For two colors the ensemble average is over SU (2)-valued gauged fields weighted 
by the usual gluonic action. The Dirac operator for quark mass m and chemical potential 
(jl is given by 

D — j u D u + m + [ijo- (2.1) 

Here, D v = d v + A v is the covariant derivative, 7^ are the Euclidean Dirac matrices, and 
A v are the SU (2)-valued gauge fields. Many of the differences between QCD with gauge 
group SU(2) and fundamental quarks and QCD with three or more colors and fundamental 
quarks originates from the pseudo-reality of SU(2). It manifests itself through the anti- 
unitary symmetry of the Dirac operator D (C is the charge conjugation operator and r 2 
acts in color space) 

Dr 2 C l5 = r 2 C l5 D*. (2.2) 

The reality of the fermion determinant is a consequence of this symmetry. However 
in our context a more important consequence is that for Nf fundamental quarks the 
flavor symmetry group is enlarged to SU(2Nf) which is sometimes referred to as the 



Pauli-Gursey symmetry [24]]. In a basis given by left-handed quarks, qz, and conjugate 
right-handed anti-quarks, a 2 r 2 q* R) (<r 2 acts in flavor space) this enlarged symmetry group 
acts on the flavor indices of the spinors 



* = ( qL „ J as * -»• Vtf, with V G SU(2N f ). (2.3) 

This symmetry becomes manifest by writing the Lagrangian for QCD with two colors in 
terms of these spinors: 

Lqcd = i^a u (D u -fiB5 0v )^- (^ T a 2 r 2 M^ + h.c.) , (2.4) 
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where cr„ = (— i, a^). The baryon charge matrix denoted by B is given by 

(2.5) 

and the matrix that includes the mass term and the diquark source term is defined by 

M = \jm 2 + j 2 (M cos + Jsin0), (2.6) 
where tan0 = j/m and 

M . (_» J), a nd °) -M. (2.7) 



Here and in the rest of the paper we consider only even Nf. The matrices M and J 
correspond to the mass term and to the diquark source, respectively. It is clear that only 
the kinetic term has the full SU(2Nf) symmetry whereas it is broken explicitly by the 
chemical potential and the source terms. 

What is more important is the spontaneous breaking of chiral symmetry by the for- 
mation of a chiral condensate. It breaks the flavor symmetry in the same way as the mass 
terms, i.e. according to SU(2Nf) — > Sp(2Nf). Along with the chiral condensate it is 
possible to form a diquark condensate. This condensate preserves the same symmetries 
as its source term, i.e. Sp(Nf) x Sp(Nf). At nonzero chemical potential and mass the 
symmetry is SU(Nf) x Ub{1)- The diquark condensate breaks this symmetry sponta- 
neously according to SU(Nf) x Z7b(1) "~ * Sp(Nf). The baryon number is no longer a 
conserved quantity. 

Although the source terms and the chemical potential break the flavor symmetry, the 
full flavor symmetry can be retained if the source terms are transformed according to 

M -> V*MV\ B -> VBV ] . (2.8) 

This global SU(2Nf) symmetry of the Lagrangian can be extended to a local flavor 
symmetry by introducing the vector field B v = (B, 0, 0, 0) with transformation properties 

iS0 

B v -> VB V V ] - -Vd u Vl (2.9) 
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Because of the spontaneous breaking of chiral symmetry, the low-energy limit of QCD 
with two colors is a theory of weakly interacting Goldstone bosons. The Goldstone man- 
ifold is given by SU(2Nf)/Sp(2Nf). If the chiral condensate is denoted by S it can be 
parameterized by 

s = um T } (2.10) 

where 



U = exp(ffl/2F) and IT = n a X a /^2N f . (2.11) 

Here, F is the pion decay constant, the fields 7r a are the Goldstone modes and the X a 
are the 2Nj — Nf — 1 generators of the coset SU{2Nf)/ Sp(2Nf). They obey the relation 
X a S = SXj and are normalized according to TrX a Xt, = 2Nf5 a b. 

The Lagrangian of the Goldstone modes is obtained by the requirement that it has 
the same symmetry properties as the underlying microscopic theory. The Lagrangian for 
E should be invariant under the local transformations 

S -> VEV T } V G SU(2N f ) (2.12) 



if the source terms and B u are transformed according to ( |2.8| , |2.9| ). We adopt the usual 
power counting of chiral perturbation theory that the square root of the quark mass, the 
square root of the diquark source, the chemical potential and the momenta are of the 
same order. Then to order p 2 it is simple to write down a Lagrangian that is invariant 
under the global transformations fl2.8| ). Invariance under the local transformation Q2.9| ) 
can be achieved by introducing the covariant derivative Jl| 

V„Tt = djtf + /i(5jSt + YtB v ). (2.13) 

The effective Lagrangian to leading order in the momentum expansion having the required 
invariance properties is thus given by 

£ (2) = ^! T r [v.SV.St - x f S - X S f ] , (2.14) 
where we have introduced the source term 

X=§^MK (2.15) 
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To leading order in the chiral expansion G is related to the quark- ant iquark condensate 
at j = fi = T = according to (ipip)o = 2NfG, and the pion mass is given by the 
Gell-Mann-Oakes-Renner relation: M 2 = mG/F 2 . 

For zero diquark source, the leading-order Lagrangian ( |2.14j) describes two different 
phases: a normal phase and a phase of condensed diquarks. They are separated by a 
second order phase transition at \i c = M/2. The chiral condensate in both phases can be 
parameterized as |2| 

S(a;) = cosaS c + sinaS d , (2-16) 

where 



a = if \i < \i c 

COS a = if yU > yU c 



;2.17) 



and 




I and S d = " ; r . (2.18) 



A non-zero value of a corresponds to diquark condensation . The orientation of the con- 
densate ( 2.16Q minimizes the static part of the leading-order effective Lagrangian (|2.14j ). 



At nonzero a it is sometimes useful to introduce rotated generators of SU(2Nf) / Sp(2Nf) 
defined by 

X a (a) = V a X a (a = 0)Vj with V a = exp(-|s d S c ). (2.19) 



At nonzero j the diquark condensate is always non- vanishing. In this case, the saddle 
point is still given by S = cosaS c + sinaE^, but a is determined by the saddle point 
equation 

4/i 2 cos a sin a = M 2 sm(a — 0), (2.20) 

where M 2 = Gy/m 2 + j 2 /F 2 . Also in this case it may be useful to introduce rotated 
generators. 

The spectrum contains 2N 2 — Nf — 1 pseudo-Goldstone bosons. In the normal phase, 
a = 0, they can be distinguished by their baryon charge. There are N 2 — 1 "usual" pions 
with charge zero, denoted by P, Nf(Nf — l)/2 diquarks with charge +2, denoted by Q, 
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and Nf(Nf — l)/2 antidiquarks with charge —2, denoted by Q>. The inverse of their 
propagator is respectively given by @, 



D p 
D Q 



p 2 + M\ 

( p 2 + M 2 _ 4/i 2 



4i/ip 
p 2 + M 2 — 4/i 2 



(2.21; 



In the phase of condensed diquarks the pseudo-Goldstone modes of the normal phase 
are mixed into four different types. We now have Nf(Nf + l)/2 P^-modes, (N 2 — Nf — 2)/2 
Pa- modes and Nf(Nf — 1) Q-modes. The inverse of their respective propagators reads 



D Ps 
D Pa 
D Q 



p 2 + M 2 + ~Ml 



P 2 + Ml + \Ml 



p 2 + M 2 iM 3 p 
iM 3 p p 2 + M 2 



[2.22) 



where 



Ml = M 2 cos(a - 0) - A/j, 2 cos 2a 

M 2 = M 2 cos(a - <p) - 4/i 2 cos 2 a (2.23) 
M 2 = 16/x 2 cos 2 a. 



By diagonalizing the inverse propagator for the Q-modes we find Nf(Nf — l)/2 Q-modes 
and Nf(Nf — l)/2 Q^-modes. The Q-modes are the true massless Goldstone modes of the 
superfluid phase at j — and fi > M/2 p. 



3 Next-to-Leading Order Effective Theory 



We wish to examine the effect of all one-loop diagrams to the free energy. Such con- 
tributions are of 0(p 4 ) in the power counting of Chiral Perturbation Theory, that is 
next-to-leading order in the momentum expansion |26|] . In this section we shall construct 
and renormalize Chiral Perturbation Theory to 0(p A ) for QCD with two colors and fun- 
damental quarks. We will closely follow the work of Gasser and Leutwyler for QCD with 
three colors and fundamental quarks |26f . 

At next-to-leading order in chiral perturbation theory, the partition function contains 
four different contributions: the tree graphs of the 0(p 2 ) effective Lagrangian & 2 \ the 
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one-loop diagrams from C^, the tree graphs of the 0(p i ) effective Lagrangian and 
the contribution from the axial anomaly ||26|| . The free energy at next-to-leading order 
can thus be written as 

S = 1S2 + Si— loop 

+ S A + S A , (3.24) 

where S 2 = J dxC^ 2 \ 5i_i 0O p is the one-loop contribution to the free energy from the 
one-loop diagrams of C^, S4 = / dxC^\ and Sa is the Wess-Zumino-Witten functional 
which reproduces the axial anomaly. 



3.1 Operators at Next-to-Leading Order 



For QCD with three colors and Nf quarks in the fundamental representation the most 
general effective Lagrangian to 0(p 4 ) has already been constructed in the literature | ]2"6] , 
27fl . The 0(p A ) effective Lagrangian for QCD with two colors and Nf quarks in the 



fundamental representation can be found along the same lines. Although the Goldstone 
manifold is different in the two cases, the same operators appear in the Lagrangian. The 
0(p 4 ) effective Lagrangian contains all the Lorentz invariant operators of 0(p A ) that are 
invariant under local SU(2Nf) flavor transformations. Keeping only the terms that are 



relevant to our case, we deduce from [26, 27| the effective Lagrangian 



C {4) = -L Ti[v u i:V T ^V^V T ^] -L^Trfv^SV^S 1 " 



-L 2 Tt 
+L 4 Tr 



Tr 



V^EV r S 



xEt + ExlTrfv.EV^t 



L 7 Tr 



-L 3 Tr 
+ L 5 Tr 



(3.25) 



-L 6 (Tr[ X £t + £ X 1 



The low-energy coupling constants L iy i = 1, . . . , 8, and H 2 are 'bare' coupling constants. 
The term related to the coupling constant H 2 is a contact term that, in our case, will 
enter only in the vacuum energy and in the quark- ant iquark condensate. Because we only 
consider quarks with equal masses, the coupling constant L 7 will not appear in any of our 
results. 



3.2 Renormalization to One Loop 

In order to compute the one-loop contributions from C^- 2 \ we expand in the fluctu- 
ations around the solution of the classical equation of motion following the analysis of 
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Gasser and Leutwyler [ 26| . The one-loop contributions from are UV divergent. They 



have to be regularized in some way, and the regularized divergences must be canceled by a 
renormalization of the coupling constants that appear in £( 4 \ If the theory is renormaliz- 
able the values of the counter-terms must be independent of M, and j, and, especially, 
the renormalization constants must be identical in the two phases. We show that this is 
indeed the case. 

The solution of the classical equation of motion in the normal phase can be written as 

E = UjlUl, (3.26) 

where U a belongs to the Goldstone manifold SU(2Nf)/ Sp(2Nf) with Sp(2Nf) subgroup 
that leaves E invariant. Using that the E can be viewed as a rotation of E c = /, i.e. 
E = V a Yi c Va with V = exp(— aE d E c /2), we can also parameterize the Goldstone manifold 
as follows 

E = V a UV^V*U T V^ = V a UlU T Vl = VJJ 2 lVl. (3.27) 

The last equality holds because U = U a=0 . The expansion around E is obtained by the 
substitution 

U -> Ue^, (3.28) 

where £ is a linear combination of the generators X a of the Goldstone manifold SU(2Nf) / Sp(2Nf) 
satisfying XJ = IXj, 

2Nj-Nf-l 

C= £ CX a , (3.29) 

a=l 

In other words, £ is a 2Nf x 2Nf hermitian traceless matrix that satisfy C.I = . The 
expansion about E can thus can be written as 

Y.J' [l + i£-^e + ---)uiVT. (3.30) 



The 0(p 2 ) effective Lagrangian is given by ( |2.14 ). Therefore, to second order in the 
fluctuations, we find 



£(2) 



— Tr[V,EV,Et - X £t _ ^t] 
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+ ^Tr[K(^t/)K(f/W) - ~V?(C/ 2 )V?(*7t£ 2 C/t) (3.31) 

- ^(c/^v^c/) + ^ 2 (t//x at t/ - tf x a iu% 

where the superscript a indicates that the original sources have been replaced by the 
rotated sources: 



A 



a = V^xV: and BS = V£B v V a , (3.32) 



and the covariant derivative V" is given by V p defined in (|2.13 )with B v replaced by B". 
Notice that B* = IB* I. 

Let us denote the classical action by S 2 - Up to second order in the fluctuations the 
action can be written as 

J dxc® = s 2 + ^-J dxTi (<u<u - [a„, e] [a„, e] - , (3.33) 

where 

r„ = ltfd v u-\ud„tf-%tfB2u-!±UBftf 

A v = ltfV u (U 2 )tf = ~UV u (tf) 2 U (3.34) 

Notice that with the above conventions V V {U£JJ) = U (d u £ + {A U ,£})U, and that the field 
strength associated with Y u reduces to Y^ u = d^Y u — d v Y ^ + \Y Y v ] = — [A M , A u ]. In terms 
of the components of £ = J2 a £, a X a where X a are the generators of SU(2Nf) / Sp(2Nf) 
with TrX a X b = 2N f S ab , we obtain 

/ dxC {2) = S 2 - N f F 2 {£, D^), (3.35) 

with scalar product defined by (/, g) = J2al dxf a (x)g a (x), and the differential operator 



(3.36) 
(3.37) 



jiven by 






d„d„C + a ab 


duC = 


d u C + Y a X 


-pab 


1 Trf[X 


V 


2N f y [ 


a ab = 


1 Trf[X a , 
2Nf v L ' 



a vbi 



(3.38) 



l-Tr^K] [X b ,A„]) +-L Tl ( a{X «,X b }). (3.39) 
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Notice that the field strength associated with f M is simply given by T ab = — |Tr([X a , X b ]T 



flU J 



The one-loop correction to the free energy is obtained by performing a Gaussian inte- 
gral resulting in 



Si-ioop = -lndet-D. 



(3.40) 



As is usual in ChPT p6f we choose to regularize the determinant by means dimensional 



regularization ||28|| . There are ultraviolet divergences that produce poles in d. The stan- 
dard calculation leads to ^] 



SWoop = E/ d d x(^5 ab 



An(d-2) 



±.ab cab , ^ f ^ yiab -pfea , ^ ±.ab ~ba\ 

° 6 + (4 7 r) 2 (rf-4)tl2 r - r - + r ° I 



In order to perform the sum over the SU(2Nf) / Sp(2Nf) generators, we need the following 
two formulae 

J2 Tr[X a A]Tr[X a B] = N f Tr[AB] - N f Tr[AIB T I) -Tr[A]Tr[B) (3.42) 

a=l 

2N 2 f -N f ~l 

E Tv[X a AX a B] = -Tr[AB] +N f Tr[AIB T I] + N f Tt[A]Tr[B]. (3.43) 

a=l 

From these two formulae we easily derive the identities 



E-pab -pba 



a.b 



2(N f - i)Tr[r^r 



4{N f - lMTr^A.A^A,] - Tr[(A^A 



(3.44) 
(3.45) 



and 



a.b 



27V / Tr[(A^) 2 ] + Tr[A M A M ] + 2Tr[A„A„]Tr[A /1 A I/ 



-Tr[a]Tr[A^] - 2iV / Tr[ t rA / ,A A 



+ 



AN] 



(JV> + 1)(JV> - 2) 



2N t 



Tr[cr 2 



(3.46) 



where we have used that X a I = IX^, -TA^J = — A M , lo T l = — cr, and IT^ U I = —T^ 
r fMU . Therefore using fl3.34|) we find 



1^^ + -^ a } 



iVf - 1. 



a, 6 



48 



-Tr 



32 



(Tr 



V^EV^Et" 
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+ 16 Tr 

-±T r 
16 

n 



V.SV.EtlTrfv.EV,^ 



+ 



2N f + 1 



48 



TV 



(3.47) 



X^ + MTrfv.SV.St 



N f 



-Tr 



+ 



^ + } ( Tr [vSt + Eyt 
32N 2 V L A A 



(7V / + l)(jV / -2) 
16JV* 



Tr 



tfvt 



/t 



The <i = 4 pole in 5i_i oop can therefore be absorbed by the following renormalization 
of the coupling constants 

Lq = Lq H — fg— A Li = L\ + jj^A L 2 = L r 2 + jgA 

L 3 = L r 3 H ^— A L 4 = L\ + -^5 = -^5 + ~8~ A 

r - T r + N f+\ j _ tt j _ T r i (JV>+i)(JV>-2) > (3.48) 



g 2 = ^+ (jVf y 2) A 



where 



A = ^A d ~ 4 + -^A d " 4 . (3.49) 

(47r) d / 2 64vr 2 v ; 

The finite counter terms in A are introduced to partially cancel the finite contributions 
of the one-loop integrals. This is the subtraction scheme traditionally used in chiral 
perturbation theory ||26|| . The cancellation of the finite counterterms will be discussed 
in detail below. In terms of the renormalized coupling constants L\ and H^, the sum 
Si-loop + ^4 remains finite for d — > 4 both in the normal phase and in superfluid phase. 
Notice that the renormalized low-energy constants depend on the renormalization scale 
A. The numerical values of these coupling constants are not yet known, but they can be 
obtained, at least in principle, from lattice QCD simulations. For QCD with three colors 
and quarks in the fundamental representation, the experimental result for the physical 
value of these coupling constants, i.e. L\ at A = m n , is that they are of the order of 10~ 3 . 
We expect that this is a reasonable estimate for their value for QCD with two colors and 
quarks in the fundamental representation as well. 



3.2.1 Wess-Zumino Term 



To account for the chiral anomaly the 0(p 4 ) Lagrangian must be supplemented by a 
Wess-Zumino term. At nonzero baryon density this term can be obtained by gauging the 
Wess-Zumino term at zero baryon density KM |30[. The latter is given by 
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where 



a = (d v T) YT Y dx v , (3.51) 

and M 5 is a five dimensional domain with space-time, M 4 , as boundary. 

At non-zero chemical potential the vector field B u = 5 U $B has only one nonzero 
component and all terms in the gauged Wess-Zumino term with more than one external 
vector field vanish. We thus find that 

IVzp, B v = 5 U , B] = /i-f- / Tr \Ba 3 ] = fi-^- [ Tr [Ba iaj a k ] e^ k ^x .(3.52) 

127T 2 JM 4 L J 127T 2 JM 4 

This term is proportional to 2i times the winding number of the field. The factor 2 
arise because fx is the quark chemical potential equal to half the baryon number chemical 
potential for QCD with two colors. The lightest excitation with nonzero winding number 
is expected to be of the order of the mass of the nucleon for QCD with three colors. Such 
terms are subleading in the low-energy limit we are considering and can be ignored. 



4 Free Energy 

In the previous section we have shown that the 1-loop contribution to the free energy 
given by to the logarithm of the determinant of the propagator matrix is renormalizable. 
In this section, we obtain an explicit expression for the free energy to one-loop order. The 
Feynman graphs that enter the free energy at this order are given in fig. 1. 




Figure 1: Feynman diagrams that enter into the free energy at next-to- leading order. The dot 
denotes the contribution from C^ 2 \ and the boxed 4 the contribution from 
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4.1 Free Energy of the Normal Phase 



In the normal phase the one-loop free energy is readily found to be given by 

= -2N f F 2 M 2 - l -{N}-l)^-\N f {N f -l)^ (4.53) 
-2N f M A (8N f L 6 + 2L 8 + H 2 ), 

where M 2 = mG/F 2 is the leading order pion mass. The terms A^ and Aq represent the 
one-loop graphs with a P-mode or both Q-modes respectively. As we shall see shortly, the 
coupling constants Lj derived in previous section Q3.48Q contain the counter-terms that 
regularize these divergent one-loop integrals. 

The one-loop diagrams with a P-mode, or both Q-modes are respectively given by 

A? = -I^'<^ + M2 -^ 2 + 1 ^) = J^r(- d -)M'. (4.55) 

For the Q-modes, we simplified the integral using that 

(p 2 + M 2 — 4/x 2 ) 2 + 16fi 2 p 2 = ((p - 2*/i) 2 + f + M 2 ) ((po + 2 if i) 2 + f + M 2 ). 

(4.56) 



Notice that the divergent part of the one-loop diagrams is exactly canceled by the renor- 
malized coupling constants ( |3.4S| ). The free energy of the normal phase reads 



-2N f F 2 M 2 1 



M 2 

P 5 " 



8N f L r e + 2L\ + H r 2 + 



2N 2 -N 



f 



256n 2 N 



T -2 In 



/ 



M 2 
A 2 



(4.57) 



As expected at zero temperature, the free energy of the normal phase does not depend 
on the chemical potential. Therefore the baryon density is zero in the normal phase. The 
free energy does not depend on the renormalization scale A: The logarithmic dependence 
of the renormalized coupling constants in A cancels the explicit logarithm in (|4.57|) . Using 
the free energy ( 4.57| ) the chiral condensate in the normal phase at next-to- leading order 
in chiral perturbation theory is given by 



dm 



2N f G 1 + 2 



M 2 

P 2 " 



IN 



La + 2LI + H r 



f^6 



2N 2 



N 



f 



128tt 2 A 



/ 



1 M 2 
- ln A^ 



(4.56 
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4.2 Free Energy of the Phase of Condensed Diquarks 



In the phase of condensed diquarks, the free energy is found to be given by 



= -2N f F\-(Mi 



Mj) 



-Ml 
4 3 



1) - 2)A P - - -N f (N, - 1)A? 



2N f (Mf - M 2 ) 2 (L + 2N f L 1 + 2N f L 2 + P 



-AN f (Ml - Mf)(Mf + ^Ml)(2N f L 4 + L, 
-8N F (Mi + jMlf{2N f L & + L 8 ) - 2iV / M 4 ( - 2L 8 + # s 



(4.59) 



where M 2 = G^/m? + j 2 /F 2 and tan0 = j'/m. 



The one-loop diagrams that contribute to the free energy contain a Ps, a Pa or the 
mixed Q modes. Their contributions are given by Q2.22 ) 



A Ps 



A Pa 
^0 



A? 



d d p 
(2%) d 

d d p 
(2n) d 

id 



In (p 2 + M 2 + -M. 
In (p° 



r(-|)(M 2 + 



Mo 



4 



" (4vr) d / 2 

= (4^ 



T(-^)(M 2 + iM 3 2 



d/2 



d/2 



= 2 



^ln(( P 2 + M 2 )(p 2 + M 2 )+ P 2 M 

T(-i)M d - 4 d 



(M 2 - M 2 ) 2 
32^ 



Mo 



i-m 

2 



-Mj 
4 3 



d/2 



+ 



(4ir) d /< 



-{Ml -Mi 



If 



Ml + Mf + yjM? + M| + M 



2M 



M 2 + Mf y/M 2 + Mf ^M 2 + Mf + M 
M| Mf 



+0((M 2 - M 2 2 ) 4 ) + 0((M 4 - M 2 4 ) 2 ) + 0(d - 4), 



(4.60) 



respectively. The masses Mx^z have been defined in ( |2.23|) . A derivation of the power 
series expansion in Mf — Mf of A^ is given in appendix B. The divergent term can 
be computed exactly, but we succeeded only to compute the finite corrections near the 
leading-order saddle point for j -C m. One easily shows that the sum of the divergent 
one-loop terms given by 



4( 4vr )d/2 



N 2 f (2M* + 2M 4 + -M 4 + M 3 2 (M 2 + M 2 2 )) - 2(N f + 1)(M 2 2 + -M 3 2 ) 2 
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is exactly canceled by the divergent contributions of the counterterms. 



In addition to the finite terms already given in the expression for A^ , there are several 
more finite terms that arise from the d — > 4 limit of the one-loop diagrams. The first type 
of contributions we consider are of the form 



A ' " wfk Md ~ A{l + s)d/2 > 



(4.61) 



and are present in A^ s , A^ A and A®. To extract the finite terms contributing to A/ we 
express A/ as 



A, 



+0(d-4) 

-A ( T \ (1 + S) 2 + — M d ~ 4 (l - S 2 ) 2 + 0(S 3 ) + 0(d - 4). (4.62) 



The second type of one-loop contributions we have to consider in the limit d 
the form 



4 are of 



A„ = ^-Jl-M d ~ 4 (M 2 - M 2 2 ) 2 . 



( 47r )d/2 



(4.63) 



It occur only in Aq . By expanding in powers of d — 4 we find that 



A 



ii 



(4vr)<V 2 64vr 2 



M d -\M( - M£Y + 0(d - 4) 



"2 A U 



(M{ — M£) + 0(d — 4). 



(4.64) 



Thus, in our subtraction scheme, there are no finite terms resulting from such contribu- 
tions. 



Using Q4.62D one easily finds the following finite contributions from the one-loop dia- 
grams to second order in M 2 — M 2 and Mf + Mf/4 — M 2 , 



2Nj - N f - 1 Nj + N 

1287T 2 



1287T 2 



+ 



{M{ -M 2 y + 



2Nf -Nf-1 



64vr 



, (M 2 + -Mj - M 



2\2 



N 2 



+ ^- 2 (Mf-M 2 )(M 2 + l Mi-M 2 ). 



(4.65) 



The second term also includes the finite contribution to Aq given in ( 4.60Q . 
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5 Landau Ginzburg Model 



In this section we study the phase transition that separates the normal phase from the 
phase of condensed diquarks. At the mean field level and for zero diquark source, it was 
found that there is a second order phase transition at \i — M/2, where M is the pion mass 
to leading order in chiral perturbation theory. The order parameter of this phase transition 
is the diquark condensate. To leading order it was found that (ipip) = 2NfG sin a, with 



Because of the contributions of the Q-integrals, the next-to-leading order analytical 
expressions for the free energy are quite complicated. However, the neighborhood of the 
critical point can be easily studied by means of a Landau-Ginzburg theory for the order 
parameter. From the leading order result we find that a is a suitable order parameter near 
the critical point. This is quite natural since a is the rotation angle of the condensate. 
Therefore if we expand the free energy in powers of a for a < 1, we should obtain 
a Landau-Ginzburg model describing the phase transition. Notice, however, that this 
Landau-Ginzburg model results from an exact calculation within the effective theory. 
Furthermore in this approach a has to be an independent variable so that the saddle 
point equation cannot be used to express a in terms of the pion mass and the chemical 
potential. We first illustrate the usefulness of the expansion of the free energy in powers 
of a by analyzing the theory at leading order and then consider the free energy at next- 
to-leading order. 

At the mean field level, the free energy in the diquark condensation phase just above 
the critical chemical potential, i.e. for a <C 1, is given by 



where the diquark source enters through <ft = tan( j/m) ~ j/m for j <C m. This free 
energy can be studied as a Landau-Ginzburg model. At <fi = there is a second-order 
phase transition where the coefficient of the a 2 -term vanishes, that is at fi = M/2. For 
zero diquark source, the order parameter obtained from the minimum of the effective 
potential is given by 



At the critical point for nonzero diquark source we find the usual mean field value of the 
critical exponent with a given by 



a = arccos(M 2 /V) @. 



n = -N f F 2 {2M 2 + 2M 2 0a + (4/i 2 - M 2 )a 2 



1 



(16/i 2 - M 2 )a 4 ) +0(a 6 ), (5.66) 
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(5.67) 




(5.68) 
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Notice that this relation implies that <f> <gi a <gi 1. Therefore the diquark condensate near 
the critical point is given by 



on 



i on 



dj i=o m d(j) 



(W,) = -— _ = _ = 2N f Ga (5.69) 



<t>=o 



where we have used that to leading order mG = M 2 F 2 . To leading order in a (// — M/2) 
expansion this result agrees with the exact mean field theory result. 

At next-to-leading order the situation is more complicated. We study the free energy 
for a chemical potential near the leading-order critical point and define p by /i — M/2 + 
Mp. The leading order results provides us with power counting rules for a, p and <f> near 
the critical chemical potential: if a ~ e, then p ~ e 2 and ~ e 3 . To obtain an expansion 
to fourth order in e, we use 

M? - M 2 4u 2 sin 2 a 2 1 4 „ 2 

M 2 + ±M 2 JX 1 2 1 

4 - 1 = cos a - 4>) - 1 = --a 2 + —a 4 + 0a + • • • . 
M 2 2 24 



(5.70) 



To this order the free energy is given by 



n M 2 ( M 2 \ 1 M 2 

-2N f F 2 M 2 = 1 + ao F^ + l v 1 + al F^J^ + a2 F^ , 



2 



+ ( 2 + J + (~ + aJ^A a\ (5.71) 



where the coefficients of the next to leading order corrections are given by 



a = 



2N 2 -Nf-1 M 2 
W,L l + 2Ll + Hi + ' ' (l-21n^), 



2N 2 f -Nf-l M 2 
ai = 8(2^L S + L5)- ^ m_, 

a 2 = 2(2N f Ll + Ll-AN f Ll-2Ll)--^±^-\n^ : 

N f M 2 

a 3 = 8(2^^) In— , 

a 4 = -|(2n / l; + l5) + |(2n / l5 + l5) + (l5 + 2n / l; + 2n / l5 + l5) 

JV/ — 1 Nf + 1 M 2 
H 7 „. T In ■ 



512 7 r 2 iV / 38Air 2 N f A 2 ' 

(5.72) 
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This free energy is typical of a system exhibiting a second order phase transition. 



For zero diquark source, the phase transition occurs when the coefficient of the en- 
teral in the free energy vanishes. Therefore we find that there is a second order phase 
transition at 



1 M 2 

— -02- 



2 
M 2 
'F 2 



F 2 

-^ f Ll-L l + 4N f L l + 2L l + ^ln^ 



M 4 
~F± 



(5.73) 



We compare this result to the pion mass at zero chemical potential at next to leading 
order in chiral perturbation theory (see Appendix A) , 



m 2 = M 2 1 + 



M 2 

T 2 



Nt + l M 2 
i(-2N f L r 4 - LI + AN f Ll + 2U,) +' J 



In 



64n 2 N f A 2 



(5.74) 



Therefore, the critical chemical potential is given by 



fi e = M/2 + Mfl c = m n /2. 



(5.75) 



For a nonzero diquark source we find that at the critical point 



a ' 



M 2N 

1 + [ai + 8a 4 ] — 



20 1 + 



M 2 
J 2 



40 



(2N f L r 4 + LI 



+ _(2JV>LS + L5 



+ 8(L r + 2N f L\ + 2N f U 2 + L\ 



3 

N, 



(5.76) 



6N] - 7N f 



Q4n 2 N 



f 



-In — 

192^/ ~ n I 2 



There are no logarithmic corrections in <p or a so that the critical exponent remains 
at its mean field value. This is somewhat surprising. For zero diquark source, half of 
the Q— modes are true massless Goldstone modes. They result from the breaking of 
the remaining symmetries of the microscopic theory at nonzero chemical potential and 
nonzero quark mass by the diquark condensate 0. For nonzero diquark source, these 
modes become massive. They are pseudo- Goldstone modes, with a square mass equal 
to M 2 0(3a + 80) /32 at leading order at /i = M/2 ( ggg ). Usually loop-graphs with 
pseudo-Goldstone modes produce logarithmic terms of the type found for instance in the 
free energy of the normal phase ( [4.571 ). In the superfluid phase, we would have therefore 
expected one-loop corrections to the free energy of the form a01n(0). These terms would 
produce a deviation from mean- field critical indices in (|5.76|) . However, we do not find 
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any logarithmic term neither in nor in a. This is a non-trivial result, and one should 
expect that to one-loop order the critical exponents of the microscopic theory are given 
by mean field theory. 

For zero diquark source, the value of the order parameter a is given by the minimum 
of the free energy. It reads 

2 ( M 2 , \( M 2 \ 

a 2 = 4 I 2/i+ — (a 2 + a 3/ u) I I 1 + 8-^a 4 J 



,2 



8 / m 7T \ ( m: 



~ (2N F Ll + LI) + ^(2N f L r 6 + L r 8 ) (5.77) 



N f - 1 12N 2 - llN f - 11 M 2 
+8(LS + 2N S V 1 + 2N,Ll + I& - ^ In ^ 

Therefore we find that the order parameter a vanishes for /i < m n /2 and that it increases 
continuously for higher chemical potential. The phase transition is thus second order. 



6 Number Density and Condensates 

The phase structure of QCD with two colors at nonzero baryon chemical potential can 
be characterized by three important observables: the quark number density, the chiral 
condensate and the diquark condensate. We compute them at zero diquark source. The 
number density is obtained from the free energy as follows 



n B 



dn 

djj, 



i dn 

M~djl 



= W f MF 2 (l + \a^\a 2 



= 4NfMF 2 a 2 (l + 
= 4NjM 2 F 2 a 2 . 



M 2 

T 2 



4(2N f U 4 + LI) - 



N f , M 2 
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32tt 2 A 2 



(6.78) 



The chiral condensate is given by 
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Finally the diquark condensate is given by 



W> = —r 



dVL 



3=0 



1 dtt 

m d(f> 
M 2 ^ 



2N f G^ + ai —ja 

S(2N f L r 6 + LI) 



( M 2 



2N 2 - Nf - 1 M2" 
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(6.80) 



The baryon number density and the diquark condensate vanish at /i = m w /2. They 
increase continuously for /x > m n /2, whereas the quark-antiquark condensate diminishes. 
At leading order it was found that in the superfluid phase (ipip) 2 + (f^) 2 = (V^) 2 ,, where 
{ipip)o is the quark-antiquark condensate in the normal phase. This relation does not hold 
at next-to- leading order. Other than that, the qualitative behavior of these observables 
is the same at leading order and at next-to-leading order. 



7 Conclusions 



As is the case for QCD with three colors, the low-energy sector of QCD with two colors 
and quarks in the fundamental representation is a theory of weakly interacting Goldstone 
bosons. What distinguishes QCD with two colors from QCD with three colors is that we 
have both mesonic and baryonic Goldstone bosons (also known as diquarks). Therefore, 
the phase transition to a phase of condensed diquarks can be described entirely by means 
of a low-energy effective theory which is completely determined by the symmetries of 
the microscopic theory. In earlier work, in which this theory was analyzed at the mean 
field level, it was found that such phase transition takes place at a baryonic chemical 
potential equal to half the mass of the Goldstone bosons. Above this transition point 
the chiral condensates rotates into a diquark condensate with a rotation angle that is 
determined by the chemical potential. Meanwhile, this result has been observed in lattice 
QCD simulations. 
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In this article we have analyzed the one-loop corrections of the leading order effective 
Lagrangian as well as the tree-graph contributions of the next-to-leading order terms in the 
effective Lagrangian. We have shown that the theory is renormalizable, i.e. all infinities 
generated by the one-loop diagrams can be absorbed into a background field independent 
redefinition of the coupling constants of the next-to-leading order terms in the effective 
Lagrangian. We have derived the counter-terms for a general background field, and by 
an explicit calculation, we have verified that they cancel the one-loop divergences both in 
the normal phase and in the phase of condensed diquarks. 

The one-loop corrections of the effective theory do not qualitatively change the predic- 
tions obtained from the mean field analysis. In particular we have obtained the physically 
satisfying result that a second order phase transition to a phase of condensed diquarks 
takes place at a chemical potential equal to half the one-loop renormalized pion mass. 
From the effective theory, we have derived a Landau- Ginsburg theory for this phase 
transition with the rotation angle of the chiral condensate as order parameter. The 
next-to-leading order corrections do not qualitatively change the mean field coefficients. 
The normal phase is characterized by a nonzero quark-antiquark condensate, a zero di- 
quark condensate and a zero baryon number density. A second order phase transition 
at n — m n /2 separate the normal phase from the diquark superfluid phase, where the 
quark-antiquark condensate, the diquark condensate and the baryon number density are 
nonzero. However, they do affect the magnitude of chiral condensate and the diquark 
condensate in a different way so that they are no longer related by the tangent of the 
rotation angle. 

At the phase transition point, half of the baryonic mesons become massless. The 
Goldstone bosons associated with the formation of a diquark condensate remain massless 
above the phase transition point. At the one-loop level these massless modes do not 
produce any infrared singularities. Neither do we find any chiral logarithms related the 
modes that become massive above the critical chemical potential. A nonzero diquark 
source does not lead to chiral logarithms either. These results corroborate with the 
absence of one-loop corrections to the mean field critical exponents. The only nonanalytic 
behavior in the free energy that we find is of the form a A ^/a 2 — 2p, (with a the rotation 
angle of the chiral condensate and p, the deviation from the critical chemical potential). 
This statement should also hold at the two-loop level. The analytic structure of the 
propagator near the critical point suggests that critical exponents remain at their mean 
field value to all (finite) orders in chiral perturbation theory. 

Our results have been derived for QCD with two colors and quark in the fundamental 
representation. Very similar low-energy effective theories have been derived for QCD 
with quarks in the adjoint representation and for QCD with three or more colors and 
fundamental quarks but with the fermion determinant replaced by its absolute value. 
The latter theory can be interpreted as QCD at nonzero isospin chemical potential and 
can be generalized to a theory with a chemical potential for each quark flavor again 
resulting in a similar low-energy effective theory. A closely related effective theory has 
been derived in condensed matter physics for disordered systems with an imaginary vector 
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potential. We expect that our conclusions for QCD with two colors will also be valid for 
each of these theories. In particular we expect the next-to-leading order corrections do 
not qualitatively change the mean field behavior. 

All our results have been derived for zero temperature. As should be the case for a 
renormalizable theory, no additional ultraviolet divergencies appear at nonzero tempera- 
ture. Just above the critical chemical potential we expect a second order phase transition 
to the normal phase at a critical temperature that is in the domain of validity of our chiral 
Lagrangian. We thus are able to obtain exact results for the critical temperature within 
the framework of chiral perturbation theory. These results as well as other properties of 
the diquark phase at nonzero temperature will be discussed in a forthcoming publication. 
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Appendix A: One-loop expressions for and F n 

In QCD with three colors and quarks in the fundamental representation, the usual 
way to compute the pion mass and the pion decay constant is to extract them from the 
axial two-point correlation function. 

For the pion mass, the Feynman diagrams of fig. 2 have to be evaluated:. 




Figure 2: Feynman diagrams that enter into the axial- vector two-point correlation function 
that contribute to the pion mass at next-to-leading order. The dot denotes a vertex from C^ 2 \ 
and the boxed 4 a vertex from . 



The pion mass is given by the pole of this two-point function in momentum space. How- 
ever, as will be explained next, we do not need an explicit calculation of this two-point 
function. The key observation is that the only finite one-loop contributions originate from 
the integral in the second diagram. At zero chemical potential, all the Goldstone modes 
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have the same mass. Therefore the one-loop diagram is proportional to 



d d p 



(2n) d p 2 + M 2 



2 (4ir) J / 2 



d-4 



(A.81) 



Because our theory is renormalizable, in our subtraction scheme divergent contributions 
of the type given in the last line are cancelled by contributions from the counter terms. 
The pion mass is thus completely determined by the tree graphs of and C^. 

The correction to the pion mass from the tree-graphs is given by twice the co- 
efficient of the quadratic term in the pion fields in . It is essential to impose the 
Euclidean on-shell condition p 2 = —M 2 in the 0(p 4 ) terms of the effective Lagrangian 



[p6[ . For example, 



Tr[d v Ed„tf] = -MVvr a + . . . 



(A.82) 



The on-shell condition results in a contribution from the L4 and L5 terms in even- 
though [L — 0. Calculating all traces in CS^ to order II 2 we find 
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(A.83) 



The pion mass is therefore given by 



4M 
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I = M I 1 + — [-2N f L 4 - L 5 + AN f L 6 + 2L 8 ] + tadpoles. 



(A.84) 



As argued before, our subtraction scheme is such that the tadpoles and the divergent 
counter terms in the combination (|3.49| ) cancel. We thus find the renormalized pion mass 
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(A.85) 



with renormalized coupling constants defined in ( |3.48| ). 

We can proceed in the same way for the pion decay constant. The contributions from 
the tree diagrams of and are easy to obtain. They are given by the coefficient 
of the d fl 7i a d ti 7r a -term in the effective Lagrangian. The contribution from the one-loop 
diagrams are again proportional to ( A. 81 ). We can therefore apply the same method as 
the one we used to determine the pion mass. We find that the pion decay constant at 
next-to-leading order reads 
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(A.86) 
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Appendix B: Q-Modes in the Diquark Condensation Phase 

In this appendix we evaluate the one-loop contribution of the Q-modes to the free 
energy. We thus analyze the following integral 



A? = -J^ln{(p 2 + Ml)tf + Ml)+plMi), 



where M\ , Mf , and M| are defined in (|2.23| ). The pole part of this integral at d = 4 can 



be easily obtained by expanding the logarithm in inverse powers of the momentum, 
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1 (M 2 + Mlf 1 p 4 M 3 4 p 2 Ml(M? + Mf 



2 p 4 2 p s p 6 

+ ••• (B.87) 

The higher orders in the expansion in inverse powers of the momentum vanish in dimen- 
sional regularization. The pole term in d — 4 is thus given by the integral 



PP(A$) = -¥±± r p d - 5 dp r d9sm d - 2 9 
(27r) d Jo Jo 



iJa P 



MfM 2 - -(Mf + M£) A - M 2 (M{ + M|) cos 2 6 - -M 3 4 cos 4 9 

(B.88) 

where we have introduced polar coordinates with po — pcos6, and the area of the d- 
dimensional unit sphere given by 

2 7r ( a! - 1 )/2 

= y^v- (B.89) 
The integrals in ( |B.88| ) can be easily calculated using the formulae 
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For the j9-integral one can derive the limiting relation for d — > 4, 



PP (f = T(- d -). 



(B.92) 



The final result for the pole part of A^ is thus given by 



(B.93) 



The finite contributions to A^ can only be obtained close to the transition point 
(a — > 0, ji — > M/2 and j — > 0). As starting point the integral is rewritten as 



M 2 + M | 



(M 2 - M|) 



We wish to compute the Q-mode contribution to the free energy up to 0(a 6 ) and to first 
order in the diquark source j, i.e. 0. To this order the integral is given by 
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The first integral can be simplified by shifting the pcrintegration. Introducing <i-dimensional 
polar coordinates as in the calculation of the pole part we then find to order a 6 . 
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The final result is 



A? 



( 47r )rf/2 



i(M 2 + M 2 ) + iM 



; ' 2 + Ill).^M d - 4 (M 2 -M, 



(47r) d / 2 8' 



I) 2 
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2\2 



In 



Mi + Ml 



2M 



{Ml - M %) 

327T 2 

Ml + M| + M 2 + M 3 V M 1 + M 2 



M 2 + M| + M| 



(B.96) 



+0((M 1 2 - M 2 2 ) 4 ) + 0((M 4 - M 2 4 ) 2 ) + 0(d - 4). 



This result correctly reproduces the pole part of Aq which was obtained at the beginning 
of this section. 



References 

[1] J.B. Kogut, M.A. Stephanov, and D. Toublan, Phys. Lett. B 464 (1999) 183-191. 

[2] J.B. Kogut, M.A. Stephanov, D. Toublan, J.J.M. Verbaarschot, and A. Zhitnitsky, 
Nucl. Phys. B 582 (2000) 477. 

[3] E. Dagotto, F. Karsch, and A. Moreo, Phys. Lett. B 169 (1986) 421; E. Dagotto, A. 
Moreo, and U. Wolff, Phys. Rev. Lett. 57 (1986) 1292; Phys. Lett. B 186 (1987) 395. 

[4] S. Hands, J. B. Kogut, M. Lombardo and S. E. Morrison, Nucl. Phys. B 558 (1999) 
327; S. Hands and S.E. Morrison, |hep-lat/99020T2| , |hep-lat/9905021| . 



[5] S. Hands, I. Montvay, S. Morrison, M. Oevers, L. Scorzato, and J. Skullerud, Eur. 
Phys. J. C 17 (2000) 285-302; 

R. Aloisio, V. Azcoiti, G. Di Carlo, A. Galante, A.F. Grillo, |hep-lat/0007018| and 
Phys. Lett. B 493 (2000) 189-196; 



Y. Liu, O. Miyamura, A. Nakamura, and T. Takaishi, hep-lat/0009009 



S.J. Hands, J.B. Kogut, S.E. Morrison, and D.K. Sinclair, Ihep-lat/OOlUO"^ 



J. B. Kogut, D. Toublan and D. K. Sinclair, |hep-lat/01040T0l ; 

B. Alles, M. D'Elia, M. P. Lombardo and M. Pepe, Nucl. Phys. Proc. Suppl. 94 (2001) 
441; 

E. Bittner, M. Lombardo, H. Markum and R. Pullirsch, Nucl. Phys. Proc. Suppl. 94 
(2001) 445; 

B. A. Berg, E. Bittner, H. Markum, R. Pullirsch, M. P. Lombardo and T. Wettig, [hcp- 



lat/0007008| . E. Bittner, M. Lombardo, H. Markum and R. Pullirsch, |hep-ph/0009~T92 



J. B. Kogut, D. K. Sinclair, S. J. Hands and S. E. Morrison, |hep-lat/0105026| . 
[6] M.A. Halasz, J.C. Osborn and J.J.M. Verbaarschot, Phys. Rev. D56 (1997) 7059. 
[7] D. Toublan and J. J. Verbaarschot, Int. J. Mod. Phys. B 15 (2001) 1404. 
[8] D. T. Son and M. A. Stephanov, Phys. Rev. Lett. 86 (2001) 592. 



27 



[9] K. Splittorff, D. T. Son and M. A. Stephanov, Phys. Rev. D 64 (2001) 016003. 
[10] J. B. Kogut and D. Toublan, Phys. Rev. D 64 (2001) 034007. 

[11] Y.V. Fyodorov, B.A. Khoruzenko and H.-J. Sommers, Phys. Lett. A226 (1997) 46; 
Y.V. Fyodorov, M. Titov and H.-J. Sommers, Phys. Rev. E58 (1998) 1195. 

[12] K.B. Efetov, Phys. Rev. Lett. 79 (1997) 491; Phys. Rev. B56 (1996) 9630; A.V. 
Kolesnikov and K.B. Efetov, Waves in Random Media 9 (1999) 71. 



[13] G. Akemann, |hep-th/0 106053 



[14] J.C. Osborn, D. Toublan and J.J.M. Verbaarschot, Nucl. Phys. B540 (1999) 317. 

[15] P.H. Damgaard, J.C. Osborn, D. Toublan and J.J.M. Verbaarschot, Nucl. Phys. 
B547 (1999) 305. 

[16] M. E. Berbenni-Bitsch, M. Gockeler, H. Hehl, S. Meyer, P. E. Rakow, A. Schafer and 
T. Wettig, Phys. Lett. B 466 (1999) 293. 

[17] P. H. Damgaard and K. Splittorff, Nucl. Phys. B 572 (2000) 478. 

[18] K. Takahashi and S. Iida, Nucl. Phys. B573 (2000) 685. 

[19] T. Guhr and T. Wilke, Nucl. Phys. B 593 (2001) 361. 

[20] P. H. Damgaard, |hep-lat/0l05TJTU. 



[21] B. Vanderheyden and A. D. Jackson, |hcp-ph/0102064j . 

[22] A. Smilga and J. J. Verbaarschot, Phys. Rev. D 51 (1995) 829. 

[23] D. Toublan and J.J.M. Verbaarschot, Nucl. Phys. B560 (1999) 259. 

[24] W. Pauli, Nuovo Cimento 6 (1957) 205 ; Giirsey, Nuovo Cimento 7 (1958) 411. 

[25] R. F. Alvarez-Estrada and A. Gomez Nicola, Phys. Lett. B 355 (1995) 288 [Erratum- 
ibid. B 380 (1995) 491]. 

[26] J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142; Nucl. Phys. B 250 (1985) 
465. 

[27] J. Bijnens, G. Colangelo and G. Ecker, Annals Phys. 280 (2000) 100. 
[28] M.J.G. Veltman, Diagrammatica, Cambridge Lecture Notes in Physics 4 (1994). 
[29] Z. Duan, PS. Rodrigues da Silva and F. Sannino, Nucl. Phys. B 592 (2001) 371. 
[30] J.T. Lenaghan, F. Sannino, and K. Splittorff, [hep-ph/0107099 . 



28 



